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Abstract. We study the effect of separate chemical and kinetic freeze-outs to the ideal hydrodynamical 
flow in Au-fAu collisions at RHIC (^snn = 200 GeV energy). Unlike earlier studies we explore how these 
effects can be counteracted by changes in the initial state of the hydrodynamical evolution. We conclude 
that the reproduction of pion, proton and antiproton yields necessitates a chemical freeze-out temperature 
of T « 150 MeV instead of T = 160-170 MeV motivated by thermal models. Contrary to previous reports, 
this lower temperature makes it possible to reproduce the pT-spectra of hadrons if one assumes very small 
initial time, tq = 0.2 fm/c. However, the pT- differential elliptic flow, V2{pt) remains badly reproduced. 
This points to the need to include dissipative effects (viscosity) or some other refinement to the model. 

PACS. 25.75.Dw Particle and resonance production - 25.75.Ld Collective flow 



1 Introduction 

Ideal fluid hydrodynamical models have been very success- 
ful in describing the bulk behaviour of particles formed in 
heavy-ion collisions at RHIC. The low-p^ single particle 
spectra as well as the transverse momentum dependence 
of elliptic anisotropy (w2(pt)) are reproduced nicely [T, 
[2] . This success has been one of the reasons to conclude 
that partonic state of matter with exceptionally low shear 
viscosity has been formed at RHIC [3j . 

However, these results have been achieved using ideal 
fluid hydrodynamical models, which assumed chemical equi- 
librium until the very end of the evolution of the system. 
This assumption is questionable, since the cross-sections of 
inelastic, particle number changing processes, are smaller 
than the cross-sections of elastic and quasi-elastic pro- 
cesses. Thus it is natural to assume that the system would 
still maintain local kinetic equilibrium when it begins to 
deviate from local chemical equilibrium. 

This kind of approach is also supported by experimen- 
tal data. The flnal hadron abundances in Au-|-Au collisions 
at RHIC can be well described by a hadron gas in approx- 
imate chemical equilibrium at Teh ~ 160-175 MeV 
The reproduction of the slopes of particle distributions us- 
ing a blast- wave model requires much lower temperatures 
around Tkin « 90-130 MeV (depending on centrality) W . 
The hydrodynamical models assuming chemical equilib- 
rium usually require freeze-out temperatures close to the 



blast-wave fit^ and cannot reproduce all observed parti- 
cle yields simultaneously. It is thus reasonable to postulate 
two separate freeze-outs: first a chemical freeze-out where 
the yields of various particle species are fixed (frozen out) 
and somewhat later a kinetic freeze-out where the particles 
scatter for the last time and the momentum distributions 
cease to evolve. 

The formalism to describe such a chemically frozen 
hadron gas has been known for quite a long time ^ . There 
have been several applications of this formalism to the hy- 
drodynamical description of heavy ion collisions at SPS 
{^/itx = 17 GeV) and RHIC (V^aa = 130 and 200 GeV) 
energies [TU1[TT[[T^[T^ . but the results have been unsatis- 
factory. The conclusion of these studies has been that if 
the chemical composition of the hadron gas freezes out at 
hadronization, an ideal fluid hydrodynamical model can 
reproduce neither the single particle spectra nor the px- 
differential anisotropy [TUfTS] . However, these studies are 
lacking in such a sense that they used the same initial state 
for both the chemically equilibrated and chemically frozen 
description. It is known that pT-distributions are sensitive 
not only to the equation of state and kinetic freeze-out 
temperature, but also to the initial pressure gradients, i.e. 
to the initial density distribution [TH]. Here we redo the 
hydrodynamic calculations once more to explore whether 
there is such an initial state which leads to an acceptable 
reproduction of the p^-spectra and elliptic flow of pions 
and antiprotons. 



" Present address: Department of Physics, Purdue Univer- 
sity, West Lafayette, IN 47907, USA 



^ With the notable exception of refs. [313. 
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2 Equation of State 

As a baseline, we use an equation of state (EoS) in chem- 
ical equilibrium. We construct it in the usual way: The 
hadronic phase is described by ideal resonance gas con- 
sisting of all hadrons and resonances listed in the Parti- 
cle Data Book [17 with mass below 2 GeV. The plasma 
phase is described by an ideal massless parton gas consist- 
ing of gluons and three quark flavours. There is a first or- 
der phase transition between these two phases at Tc = 170 
MeV. Note that this EoS is shghtly different from the EoS 
with a first order phase transition used in ref. [H] and EoS 
Q used in ref. [19j. The phase transition temperature in 
this work is = 170 MeV instead of 165 MeV and the 
number of quark flavours in the plasma phase is 3 instead 
of 2.5. The latter results in slightly larger latent heat. As 
can be seen by comparing the results here and in ref. p3] , 
these differences have only a small effect on the particle 
spectra and anisotropics. In the following we refer to this 
EoS and the corresponding initial state as CE. 

We construct an EoS out of chemical equilibrium in the 
way outlined in ref. [S] and later applied in refs. [TTlfT^fO] . 
Below the chemical freeze-out temperature all inelastic, 
particle number changing processes have ceased, but elas- 
tic and quasielastic scatterings are still frequent enough to 
keep the system in kinetic equilibrium. The quasielastic 
scatterings lead to frequent formation and decay of res- 
onances, which means that the yields of resonances and 
their daughter particles, say p-mesons and pions, still are 
in relative chemical equilibrium. This approach is there- 
fore called partial chemical equilibrium (PCE). The pres- 
ence of resonances means that the actual particle num- 
ber of any particle species is not conserved after chemical 
freeze-out, but the effective particle number is. The effec- 
tive particle number is defined as Ni = Ni + J2j '^^^''-^j: 

where Ni is the actual number of particle species i, n^*^ is 
the number of particles i formed in the decay of resonance 
j including the branching ratios, and Nj is the number of 
resonances j. The sum is over all the resonances with life- 
times smaller than the characteristic lifetime of the sys- 
tem. We take the characteristic lifetime to be 10 fm/c 
irrespective of the actual lifetime of the system. The fol- 
lowing particles and resonances have lifetimes longer than 
10 fm/c and are thus considered stable: 

TT, K, 77, oj, p, n, r]', 0, A, S, E, yl(1520), ="(1530), J7, 

and the corresponding antiparticles. The corresponding 
effective potential of resonances is obtained in a similar 

fashion: iij ~ 'Ylii'^j^ V^ii where the sum is over all decay 
products of the resonance j. 

We treat each isospin state of each particle species in- 
dependently. This results in 35 conserved quantities and 
chemical potentials in the chemically frozen stage. Tabu- 
lating the EoS as a function of 36 quantities (35 chemical 
potentials and temperature) is unpractical, but avoidable 
by taking the advantage of the isentropic nature of ideal 
fluid hydrodynamics. Since entropy is conserved, the ra- 
tio of the effective particle number density [ni = Ni/V) 
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Fig. 1. Isentropic thermodynamic trajectories below 200 MeV 
at RHIC. PCE corresponds to chemically frozen calculation 
and CE corresponds to the center of the system in chemically 
equilibrated calculation. Above the phase transition temper- 
ature nucleon chemical potential = 3/^,. The curves ter- 
minate at the kinetic freeze-out temperature ((Tkin) = 120 
and 125 MeV). Note that the difference above the chemical 
freeze-out temperature Tch = 150 MeV is due to different ini- 
tial states, see sect. |3] For PCE and CE, s/ub = 453.5 and 
395, respectively, where ub is net baryon density. 

and entropy density, fii/s, stays constant on streamlines 
of the flow. The expansion thus traces a trajectory of con- 
stant fiil s in (T, {/ii}) space and we need to evaluate the 
EoS only on this trajectory. The trajectories in the (T, /ijv) 
plane of the chemically frozen (PCE) and equilibrated sys- 
tems (CE) at RHIC are shown in Fig. m 

Qualitatively our EoS is similar to the EoSs in refs. [TT| 
[T2l[T3] . The relation between pressure and energy density 
is almost identical in chemically equilibrated and frozen 
cases whereas at fixed energy density the chemically frozen 
system is much colder than a system in chemical equilib- 
rium. 



3 Initial conditions 

We use the same boost-invariant hydrodynamic code as 
in ref. [18]. In the case of chemical equilibrium (CE) we 
use the same initial conditions as in ref. [18j for the equa- 

Table 1. Initial time, phase transition temperature, chemical 
freeze-out temperature, kinetic freeze-out energy density and 
the corresponding average temperature on kinetic freeze-out 
surface used in chemical equilibrium (CE) and partial chemical 
equilibrium (PCE) calculations. 





CE 


PCE 


TO (fm/c) 


0.6 


0.2 


Tc (MeV) 


170 


170 


Tch (MeV) 




150 


ekin (GeV/fm^) 


0.065 


0.117 


(Tkin) (MeV) 


125 


120 
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Table 2. Measured and calculated particle yields at midrapidity in most central (0-5%) Au+Au collisions at -^snn = 200 GeV. 
CE stands for chemical equilibrium, Teh ~ 170 and 150 MeV denote cases where chemical freeze-out takes places at 170 and 150 
MeV temperatures, respectively, but the initial state of the system is similar to the case CE. PCE is the final partial chemical 
equilibrium result where chemical freeze-out takes place at Teh = 150 MeV and the initial conditions are adjusted to reproduce 
the observed pT-distributions. 



dN/dy 


n+ 


TT 


K+ 




P 


P 


PHENIX [31] 


286.4 ± 24.2 


281.8 ± 22.8 


48.9 ±5.2 


45.7 ±5.2 


18.4 ±2.6 


13.5 ± 1.8 


CE 


273 


273 


44 


42 


10 


5 


Teh = 170 MeV 


268 


268 


52 


50 


25 


20 


Teh = 150 MeV 


265 


265 


51 


49 


18 


13 


PCE 


278 


278 


53 


51 


18 


14 



tion of state with a first order phase transition: The ini- 
tial entropy density distribution is a Unear combination 
of the density of participants and binary coUisions in the 
transverse plane whereas the initial baryon density is pro- 
portional to the number of participants. The parameter 
values are fixed to reproduce the py-spectra of pions and 
net protons {p — p) in the most central collisions and the 
centrality dependence of charged hadron multiplicity at 
midrapidity. The initial time, phase transition tempera- 
ture and freeze-out criteria for cases CE and PCE are 
shown in table [T] Note that the freeze-out energy density 
is slightly lower than in the corresponding calculation of 
ref. [18] because the EoS is slightly different (see sect. [2]) • 

We find that the same initialization will not work in 
the chemically frozen case. To reproduce the slopes of the 
PT-distributions we increase transverse flow by making 
the initial pressure gradients steeper and by starting the 
hydrodynamical evolution earlier. An acceptable result is 
achieved when the initial entropy density distribution is 
assumed to be proportional to the number of binary colli- 
sions in the transverse plane and the initial time is taken 
to be To = 0.2 fm/c. 

This choice of initial time is bold since it requires hy- 
drodynamics to be applicable almost immediately after 
the formation time of the partons of the system, but is nev- 
ertheless plausible [T]f5]. First, in the pQCD±saturation 
method for calculating the particle production of ref. [8], 
the average energy and particle number densities are al- 
most equal to those of a system of fixed temperature. Thus 
there is no need for particle number changing processes to 
achieve thermalization. Second, for massless particles, the 
relation between pressure and energy density is e = 3P 
for any isotropic momentum distribution. Thus the use of 
hydrodynamics to describe the build-up of collective flow 
could be a reasonable assumption even if the momentum 
distribution differs from thermal equilibrium distribution. 
Third, this ideal gas EoS may be applicable very rapidly, 
since isotropization of the momentum distribution occurs 
much faster than thermal equilibration [20]l2Tj. Finally, 
the studies of plasma instabilities support the notion of 
fast thermalization of the system [22 • 

If the proportionality between the initial entropy den- 
sity and the number of binary collisions is independent of 
the centrality of the collision, the centrality dependence 
of the final particle multiplicity is not reproduced [23J. 



To correct this, we modify the parametrization given in 
ref. [23j by assuming an impact parameter dependent pro- 
portionality constant: 

Ks{to = 0.2 fm; h) = 0.26942 fm"*^ 6^ ± 10.9 fm"'' h 
+ 453.5 fm~^, 

and parametrize the initial entropy density in the trans- 
verse plane as 

s(s; To; b) = if,(ro; h) Ta (s±ib) Tb (s-ib) , 

where Ta is the customary nuclear thickness function [23] . 
The initial baryon number distribution is obtained in the 
same way, i.e. it is taken to be proportional to the initial 
entropy density. 

4 Particle yields and chemical freeze-out 
temperature 

When separate chemical and kinetic freeze-outs are dis- 
cussed in recent literature, chemical freeze-out is assumed 
to take place immediately after hadronization at T^h = 
160-170 MeV [TO l lTT I fn i fn] . Either the ratios of all parti- 
cle species are taken to be fixed at this temperature [1111121 
113] or only strange particle yields are supposed to freeze- 
out [TUj . The idea of chemical freeze-out at hadronization 
is conceptually attractive [24ll25] and to some extent sup- 
ported by thermal models which lead to temperatures [U 
[S] quite similar to the predicted phase transition tem- 
perature [26j . However, in the context of hydrodynamical 
model, we have found that assuming chemical freeze-out 
of all particle species at T « 170 MeV, leads to proton 
and antiproton yields that are too large when the model is 
tuned to reproduce the observed pion multiplicity. In most 
central collisions we obtain dNp/dy = 25 protons and 
dNp/dy — 20 antiprotons at midrapidity instead of the 
experimentally observed yields of dNp/dy — 18.4±2.6 and 
dNp/dy = 13.5±1.8 (see Tableland fig.[2l case Tch = 170 
MeV). This discrepancy is not surprising. Some thermal 
model fits to RHIC data have actually led to tempera- 
tures near or below 160 MeV rather than 170 MeV [27l 
[28] . Thermal models also tend to lead to p/n ratios which 
are slightly larger than the experimental ratio, although 
still within experimental errors [271I281 I29 1 I30] . If the model 
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is required to fit only K/ tt and p/tt ratios, the temperature 
is lower, T « 152 MeV [30]. 

As was argued in [TS], the evolution of the average 
transverse momentum per particle, and thus the observed 
Pt spectra, is sensitive to pion number changing processes. 
In a hydrodynamical model it is thus important to de- 
scribe the pion chemistry correctly and get the pion chemi- 
cal freeze-out temperature right even if that means a worse 
reproduction of strange baryon yields. Inspired by the re- 
sults of rcf. ^27], we take T = 150 MeV as the temper- 
ature where the pion number changing processes become 
negUgible and the effective pion number freezes out. We 
expect that the better reproduction of the pion to proton 
and the pion to antiproton ratios gives a more realistic 
description of the temperature evolution and therefore of 
the evolution of particle distributions. As can be seen in 
the case Tch = 150 MeV of Tabled this approach gives 
the expected results: pion, kaon, proton and antiproton 
yields are almost exactly reproduced. The strange baryon 
yields, on the other hand, become too small at Tch = 150 
MeV. 

The possible contribution from weak decays of strange 
particles causes a significant uncertainty in determining 
the chemical freeze-out temperature [4] . According to the 
PHENIX collaboration |311I32| . the data has been cor- 
rected for the feed-down from weak decays of yl's and 
17" 's, but whether the pion spectrum contains a contri- 
bution from decays, is not clearly explained. In these 
calculations we have assumed that there is no contribu- 
tion from any weak decays in the PHENIX data. The de- 
cay of thermally produced ifg's and S^^s would increase 
the pion yield at midrapidity by 40 pions and one pro- 
ton/antiproton. For the cases Tch = 170 and 150 MeV, the 
pion yield is thus at the experimental upper limit. The re- 
duction of initial entropy by 5-10% would bring the pion 
yield closer to the experimental value and reduce the pro- 
ton/antiproton yield by the same 5-10%. The proton and 
antiproton yields at Tch = 170 MeV temperature would 
be too large in this case too. At Tch = 150 MeV, the pro- 
ton and antiproton yields would be within experimental 
errors, but one can argue that slightly larger freeze-out 
temperature Tch ~ 155 MeV is favoured by the data. 

In the case PCE the same procedure leads to pion yield 
larger than the data, which necessitates a larger adjust- 
ment in the initial entropy. However, after the adjustment, 
the conclusion is the same than for the cases Tch = 170 
and 150 MeV: Even if the pion yield contains a contri- 
bution from decays, Tch = 170 MeV chemical freeze- 
out temperature leads to too large proton and antiproton 
yields and a temperature around Tch = 150-155 MeV is 
favoured. 

On the other hand, if the experimental correction for 
weak decays is not perfect and the spectra contains feed- 
down from strange baryons, the relative increase in the 
proton/ antiproton yields would be larger than in the pion 
yield. In such a case the favoured freeze-out temperature 
would be even smaller than what is suggested here. 
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Fig. 2. (a) Proton and antiproton and (b) pion (vr^) and kaon 
(K^) pT-spectra in most central (0-5%) Au-|-Au collisions at 
-ysNN = 200 GeV compared with hydrodynamical calculations 
using different chemical freeze-out descriptions (see the text). 
The data was taken by the PHENIX collaboration [31]. For 
clarity the antiproton and kaon spectra are scaled by a factor 



5 Transverse momentum spectra 

The transverse momentum distributions for protons, an- 
tiprotons, pions and kaons in the most central collisions 
are shown in fig. O As before [H], the calculation with 
chemical equilibrium (CE) reproduces the pion and kaon 
data well, but it reproduces only the slopes of proton 
and antiproton distributions. If we proceed as in refs. |12[ 
[T5] and replace the EoS by an EoS with partial chemi- 
cal equilibrium below Tch = 170 MeV, adjust only the 
kinetic freeze-out density]^ and keep everything else un- 
changed, the result is a disaster (Tch = 170 in the figure). 
Even at (Tkin) = 100 MeV average freeze-out temperature, 
the slopes of proton and antiproton spectra are far too 
steep. Decreasing the freeze-out temperature/density does 
not help the overall fit, since the pion spectrum becomes 
steeper with decreasing freeze-out temperature. This slightly 
counterintuitive behaviour was explained in ref. [15] . In 
boost-invariant expansion, transverse energy per unit ra- 
pidity, dE'r/d?/, decreases with increasing time. In par- 
tial chemical equilibrium particle number is conserved and 



Note that in the actual calculation we use constant energy 
density as freeze-out criterion for numerical simplicity. In the 
following we talk about freeze-out temperature instead. 
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Fig. 3. Transverse momentum spectra of charged pions, 
charged kaons and protons + antiprotons in 15% most central 
Au-I-Au collisions at ^snn = 200 GeV in the ultralow-pr re- 
gion measured by PHOBOS ^3? compared with hydrodynam- 
ical calculations. The systematic and statistical errors of the 
PHOBOS data have been added in quadrature. The PHENIX 
data [31] are an average of the data in 0-5%, 5-10% and 10-15% 
centrality bins and have only statistical errors. 



transverse energy has to be distributed among the same 
number of particles. Consequently (px) decreases and the 
slope of the pr distribution steepens. 

Another problem with Teh ~ 170 MeV temperature 
can also be seen in fig. O There are too many protons and 
antiprotons as discussed in the previous section. This can 
be cured by decreasing the chemical freeze-out tempera- 
ture to Tch = 150 MeV. The calculation labelled Teh = 150 
in fig. [2] is performed in the same way as Tch = 170. Only 
the EoS is changed, but the initial state is kept the same. 
When the kinetic freeze-out temperature is kept the same, 
(Tkin) = 100 MeV, as for Tch = 170, the particle distri- 
butions are closer to the data, but the slopes are still too 
steep. The flattening of the spectra is easy to understand. 
Between T = 170 and 150 MeV particle number is not con- 
served but the energy stored in masses of heavy particles 
is converted to kinetic energy of light particles (mostly 
pions) when the system cools. Thus (px) may increase 
even if dEx/dy decreases, and (px) is larger at the time 
of chemical freeze-out when the particle number is fixed 
and pion (px) begins its slow decrease. 

Neither the steeper initial profile given by the pure bi- 
nary collision profile nor the small initial time tq = 0.2 
fm/c is sufficient alone to create large enough flow before 
chemical freeze-out to fit the data. When these two are 
used together (PCE), the fit to the data is good. The fit to 
pion data is even better than in the chemical equilibrium 
case (CE) at low px, where PCE depicts a concave curva- 
ture typical for a finite pion chemical potential ^34J . At low 
Px the calculation still suffers from excess of protons and 
especially antiprotons as if even a larger flow velocity and 
a smaller chemical freeze-out temperature were necessary. 
We have checked that the binary collision profile and short 
initial time combined with an EoS with Tch = 170 MeV 
does not lead to satisfactory reproduction of the slopes 
either, but the extra push created in the hadronic equilib- 
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Fig. 4. Pion (tt^) and antiproton (p) pT-spectra in semi- 
central to peripheral Au-|-Au collisions at ^^saa = 200 GeV 
compared with hydrodynamical calculations (see text). The 
data was taken by the PHENIX collaboration 01]. For clarity 
the spectra at centralities 20 - 30%, 40 - 60% and 60 - 80% are 
scaled by factors 10~^, 10~^ and 10~^, respectively. 



rium stage between T = 170 and 150 MeV temperatures 
is necessary to fit the data. 

Figure [3] shows the results for the very low transverse 
momentum measured by the PHOBOS collaboration [33] 
as well as the spectra measured by the PHENIX collabo- 
ration [31] at larger pt- The model PCE works well also at 
this Px region reproducing the flat behaviour of the data. 
The spectrum of pions is slightly below the data. This was 
seen already in table [2] where pion multiplicity was shown 
to be slightly below the data but still within the exper- 
imental error. The increase in the initial entropy of the 
system to increase pion multiplicity would, however, ne- 
cessitate even lower chemical freeze-out temperature not 
to exceed the observed kaon multiplicity. 

Pion and antiproton py-spectra at various centralities 
are shown in fig. [U The pattern is familiar from earlier 
hydrodynamical studies: The more peripheral the colli- 
sions, the narrower the px range where hydrodynamics 
can reproduce the spectra. An interesting detail is that 
antiproton spectra in semi-central collisions (20-30% and 
40-50% of total cross section) is slightly too flat. As men- 
tioned in the introduction, blast wave fits give higher ki- 
netic freeze-out temperatures in peripheral than in cen- 
tral collisions. Hydrodynamical models have traditionally 
used fixed freeze-out temperature/density for simplicity. 
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Fig. 5. Elliptic anisotropy of pions and antiprotons vs. 
transverse momentum in minimum bias Au-(-Au collisions at 
y^SNN = 200 GeV calculated using different chemical freeze- 
out descriptions (see text) and compared with the data by the 
STAR [35] and PHENIX [36] collaborations. 

but the overshooting of antiproton spectra here suggests 
that a slightly larger kinetic freeze-out density in semi- 
central than in central collisions might lead to better re- 
sults. 



6 Elliptic anisotropy 

Unfortunately the success in the previous section cannot 
be repeated for elliptic anisotropy. Figure [5] shows the px- 
differential elliptic anisotropy of pions and antiprotons in 
minimum bias Au-|-Au collisions. The trend is similar to 
that in ref. [T^]. The chemical equilibrium result (CE) re- 
produces the pion data excellently and is slightly above the 
proton data. Partial chemical equilibrium (PCE) leads to 
a slope of the pion anisotropy that is too large and the 
calculated result is clearly above the data. The proton 
anisotropy is reproduced at pt < 600 MeV, but increases 
too fast at larger momenta. 

It is worth noticing that the different pr-differential V2 
is mostly due to the different description of the hadron gas. 
As shown in ref. T2] , the same initial state leads to differ- 
ent V2 in chemically equilibrated and frozen cases. The ba- 
sic reason for this is in the different temperature evolution 
in the hadron gas. When both systems have reached the 
same density, the collective flow field and its anisotropics 
are almost identical in both cases. However, the chemi- 
cally equilibrated system is hotter and thus the random 
thermal motion is stronger. Thermal motion smears away 
the underlying anisotropy of the collective flow field, and, 



therefore, the final V2 of particles of the chemically equi- 
librated system is smaller than the V2 of particles of the 
corresponding chemically frozen systenjf]. 

Comparison of these results with the recent results ob- 
tained using viscous hydrodynamics [5711^155] is partic- 
ularly interesting. Dissipation reduces elliptic flow and it 
was claimed in ref. [39] that even the postulated minimal 
shear viscosity 7 = 3^ causes very large suppression. As 
shown here, suppression of elliptic flow is required to fit 
the data when chemical equilibrium is lost in the hadron 
gas. However, the actual size of the allowed viscous cor- 
rection is unknown since the EoS affects especially the 
proton anisotropy. At least in the chemically equilibrated 
case, the lattice QCD based EoS leads to even larger pro- 
ton V2{pt) than the EoS used here [TSj. Whether finite 
shear viscosity leads to a correct reproduction of the data 
remains to be seen as well as how large of a viscosity is 
allowed in plasma and hadronic phases [40] . 

7 Summary and discussion 

Contrary to previous reports [T2 | [T3 l [T4 l [T5] , we have shown 
that it is possible to reproduce both (nonstrange) hadron 
yields and their transverse momentum spectra using ideal 
fluid hydrodynamics with separate chemical and kinetic 
freeze-outs. The difference between our approach and the 
previous studies is that we used a different initial state for 
the chemically frozen system than for the system in chem- 
ical equilibrium, adjusted the initial time to be as small as 
To = 0.2 fm/c and changed the chemical freeze-out tem- 
perature to be below the phase transition temperature to 
Teh = 150 MeV. However, our success is incomplete. Even 
after these adjustments of the model, we were not able 
to reproduce the pr-differential elliptic flow of pions and 
protons. To describe particle yields, their pT-distributions 
and anisotropics, the description of expansion has to be 
amended, most probably by including dissipative effects: 
viscosity and diffusion. 

Nevertheless, our results demonstrate how sensitive 
the expansion dynamics during the hadronic phase is to 
hadron chemistry. Good knowledge of hadron chemistry 
is therefore essential if we wish to extract the viscosity of 
the QGP from the experimental data using viscous hydro- 
dynamics. The recent method of choice for this purpose 
has been the so called hybrid hydro-|-cascade model where 
the expansion is described using hydrodynamics until the 
system is hadronized and the hadronic phase is described 
using hadronic cascade [40 l l4T 1 l42 j l43] . These models have 
the advantage that the dissipative effects and hadronic 
chemistry are included in the model and the kinetic freeze- 
out is a result of hadronic cross sections without any free 
parameters. So far cascade models used in these hybrid 
models have been limited to two particle scatterings, but 
it has been argued that multiparticle processes are essen- 
tial for reproducing the proton and antiproton yields |44i 
[45] . It is therefore not surprising that hydro+cascade mod- 
els may have difficulties in reproducing antiproton spec- 

For a more detailed discussion see ref. [15j . 
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trum at the lowest values of transverse momentum |43j . 
On the other hand, multiparticle processes are included 
in hydrodynamical models and as argued in sect. [5l pro- 
cesses where heavy particles annihilate and form lighter 
particles, and vice versa, are important for the build-up 
of flow. For this reason we consider that a hydrodynamical 
description of the hadronic phase is still worth pursuing 
as a complementary tool to hybrid models, even if there 
are transport models where the multiparticle processes are 
included j45j . 

Another conclusion of these results is that the final 
spectra are sensitive to the details of the initial state. As 
already stated in ref. piT, theoretical constraints to the 
initial state and especially to its shape are of utmost im- 
portance if we wish to extract the properties of the QGP 
from experimental data. 

Some arguments favouring the use of hydrodynamics 
already at initial time tq = 0.2 fm/c were given in sect.|31 
We admit that in a certain sense our result is similar to 
ref. [13] , where Kolb and Rapp claimed that it is not possi- 
ble to fit the pT-spectra without pre-hydrodynamic trans- 
verse flow. Some of our arguments favouring the use of 
hydrodynamics at r = 0.2 fm/c centered on the fact that 
for the relation e = 3P to hold, it is sufficient that the 
momentum distribution is isotropic. It docs not need to 
be thermal. Thus one may say that we are using hydrody- 
namics to calculate the pre-hydrodynamic transverse flow, 
although the flow field in our case is different from the one 
used in ref. [13]. 

The low chemical freeze-out temperature Tch — 150 
MeV clearly below the phase transition temperature Tc = 
170 MeV leads to another problem. The chemical freeze- 
out temperature is almost independent of centrality ^30j . It 
was argued in ref. [25] that such an independence requires 
either chemical freeze-out at phase transition or extremely 
steep temperature dependence of the inelastic scattering 
rate. One way to solve this apparent discrepancy is to 
assume that the chemical equilibration in the hadron gas 
close to Tc is very fast indeed due to massive Hagedorn res- 
onances [46] . These resonances couple to multi-pion states 
and baryon-antibaryon pairs and thus lead to fast chemi- 
cal equilibration, but they appear only close to the critical 
temperature. 

The treatment of hadronic chemistry and choice of 
chemical freeze-out used in this paper is by no means 
the final one. Since the data seem to favour the notion 
that strange baryon yields are fixed at higher temperature 
than the yields of pions, kaons and nucleons/antinucleons, 
it would be relatively easy to use two separate chemical 
freeze-out temperatures: one for strange baryons and an- 
other for everything else. Another relatively straightfor- 
ward improvement is to redo the calculation of the chem- 
ical relaxation rate of pions [47] for the RHIC environ- 
ment and to use it to dynamically find the local chemical 
freeze-out temperature by comparing this rate to the ex- 
pansion rate. This approach is particularly appealing if the 
kinetic freeze-out surface is also found by comparing the 
scattering and expansion rates. Finally, it would be use- 
ful to calculate the particle yields dynamically during the 



evolution using actual rate equations for particle number 
changing processes. For baryons, preliminary study was 
already done in ref. [48j . but in that work mesons were 
supposed to be in chemical equilibrium and the EoS was 
independent of baryon and antibaryon densities. A proper 
calculation would require dynamical treatment of mesons 
and calculating the EoS using actual particle densities. 
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